Double criticality of the SK-model at T = 0. 
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Numerical results up to 42nd order of replica symmetry breaking (RSB) are used to predict the 
singular structure of the SK spin glass at T — 0. We confirm predominant single parameter scaling 
and derive corrections for the T — order function q(a), related to a Langevin equation with 
pseudotime 1/a. a — and a — oo are shown to be two critical points for oo-RSB, associated with 
two discrete spectra of Parisi block size ratios, attached to a continuous spectrum. Finite-RSB-size- 
scaling, associated exponents, and T = 0-energy are obtained with unprecedented accuracy. 
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The low temperature limit usually simplifies consider- 
ably the properties of magnetically ordered phases. Re- 
search in recent decades has however shown that frus- 
trated systems can have rich behaviour even at T = 0. 
Spin glasses [l[ are an extreme example from condensed 
matter, while others are a feature of computer and in- 
formation science in problems such as hard satisfiability 
and error-correcting codes. In particular, even the poten- 
tially soluble infinite-range Ising spin glass model of Sher- 
rington and Kirkpatrick [2| has left open many puzzling 
questions. Parisi devised an ansatz [3[ for the order pa- 
rameter of the SK-model, based on an infinite hierarchy 
of so-called replica symmetry breakings and related hi- 
erarchically to the distribution of overlaps of metastable 
solutions [101 ]. The determining equations for this ansatz 
have recently been rigorously proven to be exact Q , but 
its explicit solution remains elusive. Also only recently 
has the T = SK problem been recognized as a critical 
one-dimensional theory 0, Q . 

In view of the paradigmic role that the SK-model has 
played in the understanding and development of the sta- 
tistical physics of complex systems, together with the po- 
tential that further comprehension of its subtleties has for 
extensions to other more-complicated systems in many 
fields of science, especially those involving zero- (or effec- 
tively zero-)temperature replica-symmetry-breaking [l7| 
transitions, it seems important to pursue the better un- 
derstanding of T = RSB in the SK model. This letter is 
concerned with such a study and the exposure of several 
novel features, including new critical spectra, invariance 
points and quasi-dynamics. 

Parisi's order parameter is a function q(x, T) on an 
interval < x < 1, the limit of a stepwise function 
qi(T), Xi(T) determined by extremization of a free en- 
ergy It provides the hierarchical distribution of pure 
state overlaps P(q) through P(q) = dx/dq [ltj. Parisi's 
original work considered numerically an approximation 
with a small finite number of steps, but most recent stud- 
ies of the SK model have been based on self-consistent 
solutions for his later non-trivial continuous order func- 
tion, typically perturbatively in the deviation from the 
finite-temperature phase transition. Here the analysis is 
considered explicitly at T — using very accurate studies 



of a very large sequence of RSB orders. 

In the limit of zero temperature Parisi's order func- 
tion may be replaced by q(a),0 < a < oo, q(a) = 
lim-r^o qparisi{aT) Already, in recent work to 5RSB 
we observed a predominant single parameter scaling 
for q(a) and used this to predict a simple order function 
q(a) = V^Fa/ (2£) erf (£/a) = iFi(l/2, 3/2, -£ 2 /a 2 ) where 
a single correlation 'length' (in RSB-space) ^ w 1.13 char- 
acterized the non-trivial scaling dependence on a [l8| . 
Now, by means of redesigned numerical procedures and 
analytical transformations we have been able to improve 
the previous 10-dimensional extremization (5-RSB) up to 
an 84-dimensional extremization of the SK-energy (42- 
RSB). To the best of our knowledge, this is by far the 
highest order calculation of RSB, allowing for the first 
time finite 'size' scaling on a rather large one-dimensional 
'lattice' with up to 42 'sites', where each site stands for 
one RSB-order 19]. This has enabled many hitherto un- 
known interesting new features to be identified. 

Let us recall the SK Hamiltonian 7i = 2j<j JijVi&ji 
a = ±1, with Gaussian-distributed Jy independently 
distributed over all pairs of sites with mean zero and 
variance AT -1 . The T = 0- limit [11] of the free energy 
can be written as the extremization with respect to the 
Parisi order parameter plateaux qi and steps a.; of 

B.=Ei(Wi)-^Mnw eN ) o) 

l—l l—K 

where % = acts like a RSB-transfer matrix as 

tfriifV*) = f dh ie xp(-^(h i+1 hif) f (h^ 

with % = qi - q i+1 , r t = Oi/oj-i, qi = 1, h K+2 = 0. 

The selfconsistently obtained q- and a-level distribu- 
tions are shown in Fig[TJ High RSB-orders k demon- 
strate several interesting features. We remark first upon 
the emergence of spectral bands and the formation of a 
dense region of a-levels close to a w £. The spectral char- 
acter is further elaborated in Fig[2]by plotting the ratios 
r\ K ' = a\+i/a\^ on the unit interval < i/k < 1 for all 
RSB-orders K. Already by eye one sees that the Parisi 
blocksizc ratios approach characteristic limits as k — > oo, 
further justified by fitting to Pade approximants. These 
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FIG. 1: q 8 -spectra (left) are plotted versus RSB-order re; a-level spectra (right) are displayed by log(a* ) versus c log(re) 
with c « 4/3; asymptotic linear behavior of diverging log(a^ K ') and a dense regime for a« = 0(re°) are observed as re — > oo. 
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FIG. 2: Left figure shows results for Parisi block size ratios r^-, = a K-i+i/ a K-i (black dots) normalized to the unit interval 
< i/n < 1, connected by lines with period-10 alternating colors from 2-RSB (red, bottom) to 42-RSB (red, top); ao = 
co, a K +i = 0. Right figure shows extrapolations to three different spectra at oo-RSB: a discrete spectrum on line ((0, 0), (0, 1)), 
a continuous distribution of r = 1 along the line ((0, 1), (1, 1)), and a second discrete spectrum from (1, 1) to (1, 0). 
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FIG. 3: Discrete spectra for Parisi block size ratios n (ratios 
of largest a) and ft = r K _u_u (ratios of smallest a) as ob- 
tained for i — 0, ..,22 in the re — > co-limit (dots) compared 
with fit curves interpolating (discrete) analytical predictions 



demonstrate the existence of two discrete spectra on the 
lines d = {(0,0), (0,1)} and £ 2 = {(1, 0), (1, 1)}, re- 
spectively connected at the singular accumulation points 
(0,1) and (1,1) to a continuum on {(0, 1), (1, 1)}. Ex- 
plicitly, the continuation to k = oo yields (i) a discrete 

spectrum of the form fi = /i7n K ^ 00 7'^ i = 1 — + 1) 
for integer-valued i — 0,1,2,... on C\ and (ii) another 
discrete spectrum for r.; = Hindoo on line £2, well 

approximated by r< = (1 + f /((i + §)* - |) 2 ))" 1/2 [H, 
(see Fig 3) while (hi) continuations at fixed finite i/n 
(see FigJ5] right) demonstrate the continuous spectrum 
(as k — ► 00 ) . The two discrete spectra reflect the inequiv- 
alent non-analytic behaviors of the order function q(a) at 
its limiting values at a = 00 and a = 0. Fig0] shows our 
recent conjecture for the analytic continuation of the k- 
RSB stepped order function q^(a) = ^f-2-erf(^4= — ) 
with lim-re—joo w K = and ^ K — > ^ w 1.13 is confirmed as a 
good fit by the numerical solutions from k = 10 to n — 42 
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FIG. 4: (color-online) 10- to 42-RSB-data compared with 
analytical order function model q(a) = iFi (a, 7, — £ 2 /(a 2 + 
w))a/Va 2 + w for Pi = {a = §,7 = |,W = 0(1(T 6 )} 
(dashed-blue) and P 2 = {a « 0.53,7 « 1-71, w « 0.02} (or- 
ange); the insert zooms a tiny nonlinearity of q(a) around 
a w 0.4 with a maximum in q'(a); the data are best repro- 
duced for parameters P2 {q for Pi is too flat for small a). 



but on smaller a-scales a tiny correction is found, better 
resolved by the derivative q'(a) which shows a maximum 
near a ~ 0.344. This correction can be incorporated in 
an improved analytical model fit function [2l| by 



l F 1 (a K ,7 K ,-^/(a 2 + WK )). (2) 



The flow of the variational parameters as k — > 00 remains 
close to our previous proposal Q : the limits slightly mod- 
ified to a = = 0.53 ± 0.01, 7 = 7oo = 1.71 ± 0.02, 
£ = £oo = 1.16 ± 0.01, and w monotonically falling to- 
wards zero [22l | . 

We observe that the large a-expansion q°°(a) = 1 — 
ii4 + ^±^-^ + 0{a~ 6 ) w 1 - 0.42a" 2 + 0.16a" 4 + 

7 a z 7(7+1) 2a 4 1 V / 

0(a~ 6 ) of Eq|2] agrees well with Pankov's q(a) = 1 — 
0.41a -2 obtained from the standard oo-RSB differential 
equation and T — ► scaling Q . The deviation of a and 
7 from the rational values i and | of our original pro- 
posal Q leads to a small correction to a linear rise of 
q(a) at small a and a pronounced maximum of q'(a) at 
a = 0.344. This feature confirms low T results of Crisanti 
and Rizzo [13], obtained perturbatively. It is also inter- 
esting to compare our results for q(a) with the predic- 
tions of the Vannimenus-Toulouse-Parisi scaling ansatz 
qparisi{x,T) = f(x/T) < q max (T) @, H, combined with 
the Parisi- Toulouse hypothesis Q and the knowledge of 
the behaviour on the de Almeida-Thoulcss line. The com- 
parison is qualitatively quite good but not perfect, most 
noticeably deviating at low a, the VTP prediction being 
slightly lower and more curved. 

Our analysis also gives strong hints for the existence 
of invariance-points in a-space, i.e. points where q^ K \a) 
or one of its derivatives d™q( K '(a) does not vary under a 
change of RSB-order. They appear to be fixed points un- 
der the renormalization group which decimates the RSB- 
order [||. For example q^ K \a) does not change with k 



at a — a f» 0.401; for all calculated orders the changes 
q( K+1 )(a) _ g( K )( a ) are negative for a < a and positive 
for a > a; they decay towards zero in oo-RSB without 
change of sign, such that the contributions accumulate. 
As yet we have no physical explanation for these new 
critical "a- values". 

Our accurate results, O(10~ 12 ), for such a large num- 
ber of RSB orders also enables us to determine finite- 
RSB-size scaling exponents and the ground state energy 
to unprecedented accuracy. The energy of the SK-modcl 
is plotted versus RSB-order n in FigJS] Within the or- 
ders studied both the energy and the susceptibilities as a 
function of k have reached the asymptotic large- k scaling 
regime, fitting accurately the power law behaviors 



with E n 
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= -0.7631667265(6), z E 
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(x — 1) vanishing exponentially [23j. Figrj] also shows 
few-parameter analytic functions that yield good fits for 
the full range of RSB-orders. The results for and 
1 = x(T = 0) = J °° da a q'(a) are in accord with Eq.(2). 

In an earlier 5-RSB study of the T = distribution of 
local fields P(h) = iV" 1 £ 4 (<5(/i - £V J ij°j)) there ap- 
peared to be a deviation from linearity in \h\ at very 
small h 0. In the light of the extension to 40- RSB, we 
are able to show that this is an artefact of finite RSB or- 
der and the power law P(h) w 0.3\h\ for small h [H|,[III is 
confirmed. In any specific finite-RSB order, small oscil- 
lations are observed around P'(h) w 0.3 near h = xi but 
they shrink to zero as xi — *■ in the oo-RSB limit. Thus 
in the present context the pseudogap behavior P(h) ~ \h\ 
can be viewed as being generated by the decay of finite-K 
RSB gaps as n — * 00. 

One may further relate our results to a pseudo- 
dynamical field theory. The order function q(a) in Eq.(2) 
can be recognised for w = as the solution of the conflu- 
ent hypergeometric differential equation (a 4 /£ 2 )<9 2 g(a) + 
((3 - 27)a 2 /£ 2 - 2)a<9 Q g(a) + 4ag(a) = 0, or given by 
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log(g(a)) 



9(£/a) rv ' 50(a) ' ^ ' V2o^(£/a) 
U = \/%a sinh(0(a)) - ((7 - l/2)a/£ + f/a) <t> 2 {a) (4) 

provided sinh(0) is truncated at third order. In fact, how- 
ever, the numerical solution of the sinh(0)-model Q pro- 
vides as good a fit to our data as Eq. ^ and thus mimics 
similarly well the SK-model at T = 0. The relationship 
to dynamics follows from an inverse relationship between 
the Parisi 'length' x (or a as here) and the timescale to ex- 
plore the corresponding clusters of pure states in dynam- 
ics, with the characteristic timcscalcs increasing sharply 
as x is reduced [l5| ■ Taking 1 /a as a pseudo-time Eq. ([j) 
is identified as a (Langevin) pseudo-dynamical equation. 
The parameters a and 7, which model the order function 
q(a) analytically in between the singular points a = 0, 
a = 00, appear as two new coupling constants in the 
Hamiltonian. We also recall here the known relationship 
between replica-symmetry breaking and slow dynamics, 
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much in analogy with the Halpcrin-Hohcnbcrg theory [9] 
of slow-dynamical (and spatially-critical) behavior. 

In summary, we have obtained several new results 
for the paradigmic Sherrington-Kirkpatrick spin-glass at 
T = 0: (i) we have found that the Parisi order func- 
tion can be represented by a ID field theory with two 
inequivalcnt critical limits a — and a = oo, the order 
parameter q(a) showing critical behavior in both limits 
of long (a = 0) and short (a = oo) pseudotimes. The 
key issue for criticality is the presence of two discrete 
spectra for the Parisi block size ratios (at a — and at 
a = oo), attached to a continuous spectrum for finite a. 
(ii) All orders of RSB up to 42nd order have been solved 
to great numerical accuracy to enter deep into the largc- 
RSB scaling regime, enabling us to obtain new power 
laws, the SK-energy up to O(10 -10 ), and corrections to 
single parameter scaling in a hypergeometric order func- 



tion model. 

It would be interesting to extend to systems, both 
range-free and finite-range, and to include real temporal 
dynamics and quantum effects. As examples of poten- 
tial non-trivial range-free extensions one might consider 
K-SAT or LDPC code problems under dynamical algo- 
rithmic optimization using either classical or quantum 
dynamical annealing [24[ ■ The existence of true RSB in 
finite-range problems at equilibrium is still controversial 
but our study could potentially usefully extend current 
Ginzburg-Landau- Wilson functional integral theory [l(|. 
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